Observation of an energy-dependent difference in elliptic flow between particles and 

anti-particles in relativistic heavy ion collisions 
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Elliptic flow (1)2) values for identified particles at mid-rapidity in Au+Au collisions, measured by 
the STAR experiment in the Beam Energy Scan at RHIC at y'sjviv — 7.7-62.4 GeV, are presented. 
A beam-energy dependent difference of the values of V2 between particles and corresponding anti- 
particles was observed. The difference increases with decreasing beam energy and is larger for 
baryons compared to mesons. This implies that, at lower energies, particles and anti-particles are 
not consistent with the universal number-of-constituent-quark (NCQ) scaling of V2 that was observed 
at = 200 GeV. 
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PACS numbers: 25.75.Ld, 25.75.Nq 



Lattice Quantum Chromodynamics (QCD) predicts 
that at sufficiently high temperatures, T, and/or high 
baryonic chemical potentials, ^b, normal nuclear mat- 
ter will undergo a phase transition to a state of mat- 
ter where quarks and gluons are deconfined, called the 
Quark-Gluon Plasma (QGP) [T[. A Beam Energy Scan 
(BES) program Q has been carried out at the Rclativis- 
tic Heavy Ion Collider (RHIC) facility to study the QCD 
phase structure over a large range in T and yu^. Parti- 
cle production in heavy ion collisions with respect to the 
event plane (EP) can be characterized by the following 
Fourier expansion: 



dN 



d{(j) - 



cos - ^)] 



(1) 



where is the azimuthal angle of the particles, n the 
harmonic number, v'^^ the observed Fourier coefficient 
which has to be corrected for the EP resolution toget f„, 
and ^' the reconstructed EP azimuthal angle 0, The 
second harmonic coefficient is denoted as elliptic fiow, 

Elliptic flow measurements have been used to conclude 
that strongly interacting partonic matter is produced in 
Au-|-Au collisions at ^/sNN = 200 GeV and that V2 devel- 
ops in the early, partonic, stage. This conclusion is based 
in part on the observed scaling of V2 versus the trans- 
verse momentum, pt^, with the number of constituent- 
quarks (NCQ) for hadrons at intermediate pr (2 
to 5 GeV/c). Deviations from such a scaling for iden- 
tified hadron V2(j>t) at lower beam energies is thus an 
indication for the absence of a deconfined phase 

In a hydrodynamic picture, V2 arises in non-central 
heavy ion collisions due to an initial pressure gradient, 
which is directly connected to the eccentricity. This leads 
to particle emission predominantly in the direction of the 
maximum of the pressure gradient. During the expan- 
sion of the system the pressure gradient decreases, which 
means that elliptic flow primarily probes the early stage 
of a heavy ion collision. 

For Au-|-Au collisions at ^^snn = 200 GcV, a mass or- 
dering in V2(j>t) between the different particle species 
was observed at low transverse momenta [pT < 2 
GeV/c) 0,1 This behaviour can be described by 
non- viscous hydrodynamic calculations [TT|-[T6|. The rel- 
ative mass ordering can be suppressed by using the re- 
duced transverse kinetic energy (my — mo) instead of pr, 
with rriT = \/p\ + niQ and mo being the mass of the par- 
ticle. At large (m^ — toq), a splitting in V2{mT — mo) 
between baryons and mesons was observed which can- 
not be described by hydrodynamic calculations. This 
splitting can be explained, in part, by assuming that the 
particle production occurs via coalescence of constituent 
quarks [l7|. 

The V2 values for tt^, K^, Kg, p, p, 0, A, A, S~, 
S^, f2~, measured at mid-rapidity in minimum 



bias Au-|-Au collisions will be reported. The data were 
recorded by STAR, the Solcnoidal Tracker at RHIC, for 
= 7.7, 11.5, 19.6, 27, 39, and 62.4 GeV in the 
years 2010 and 2011 as part of the BES program 0. 

STAR is a multi-purpose experiment at RHIC with a 
complete azimuthal coverage. The main detectors used 
for the data analysis were the Time-Projection Cham- 
ber (TPC) [l3 for tracking and particle identification 
at pseudo-rapidities I77I < 1.0, and the Time-of-Flight 
(TOF) detector. A minimum bias trigger was defined us- 
ing a coincidence of hits in the Zero Degree Calorimeters 
(ZDC), Vertex Position Detectors (VPD), or Beam-Beam 
Counters (BBC) [l^. To suppress events from collisions 
with the beam pipe (radius 3.95 cm), an upper limit cut 
on the radial position of the reconstructed primary ver- 
tex of 2 cm was applied. In addition, the z-position of 
the vertices was limited to values less than ±70 cm. Col- 
lisions within a 0-80% centrality range of the total re- 
action cross section were selected for the analysis. The 
centrality definition is based on a comparison between 
the measured track multiplicity within \'q\ < 0.5 and a 
Glauber Monte-Carlo simulation |l9| . 

The particle identification and yield extraction for 
long-lived charged hadrons (p, p, tt*, K^) was based on a 
combination of the ionization energy loss, dE/dx, in the 
TPC, the reconstructed momentu m (p ), and the squared 
mass, m , from the TOF detector |20i. Short-lived parti- 
cles which decay within the detector acceptance such as 
0, A, A, S-, S , n-, n , and K° were identified using 
the invariant mass technique. The combinatorial back- 
ground to the weakly decaying particles like A and S was 
reduced by topological reconstruction. The remaining 
combinatorial background was fit and subtracted with 
the mixed event technique [20| . 

The event plane was reconstructed using the proce- 
dure described in Ref. In order to reduce the effects 
of non-fiow contributions arising mainly from Hanbury- 
Brown Twiss correlations and Coulomb interactions, the 
event plane angles were estimated for two sub-events sep- 
arated by an additional 77 - gap instead of using the full 
TPC event plane method j20l |. For such an "77-sub-EP" 
reconstruction, one uses only the particles from the oppo- 
site 77 hemisphere with respect to the particle of interest 
and outside of an additional ??-gap of \ri\ > 0.05. The 
non-flow contributions were studied for the six beam en- 
ergies by comparing different methods of extracting V2 
for inclusive charged hadrons [l9| . The four particle cu- 
mulant W2{4} strongly suppresses non- flow contributions. 
It has been shown that the difference between V2 (?7-sub) 
and V2{4:} is 10-20% for 19.2, 27, and 39 Ge V and de- 
creases with decreasing energy. All observed values (vf^^) 
were corrected on an event-by-event basis using the EP 
resolution [21| which was calculated by comparing the 
two 77-sub-EP angles [l9| . 

For each particle species, the cuts used for particle 
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FIG. 1. (Color online) The elliptic flow V2 of protons and anti-protons as a function of the transverse momentum, pt, for 
0-80% central Au-(-Au collisions. The lower panels show the difference in V2(pt) between the particles and anti-particles. The 
solid curves are fits with a horizontal line. The shaded areas depict the magnitude of the systematic errors. 
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FIG. 2. (Color online) The difference in V2 between particles 
{X) and their corresponding anti-particles (X) (see legend) as 
a function of y^sjviv for 0-80% central Au-|-Au collisions. The 
dashed lines in the plot are fits with a power-law function. 
The error bars depict the combined statistical and systematic 
errors. 



identification and background suppression were varied to 
estimate the systematic uncertainties. The errors were 
also estimated by varying the methods used to flatten 
the EP, to obtain the yields, and to extract the V2 val- 
ues. A more detailed description of the detector setup 
and the analysis can be found in Ref. [1^ . 

In Fig. [1] the pT dependence of the proton and anti- 
proton V2 is shown for Au-|-Au collisions at = 
7.7, 11.5, 27, and 39 GeV. At all energies, the V2 val- 
ues increase with increasing pT- At px = 2 GeV/c, the 
magnitude of V2 for protons increases with energy from 
about 0.10 at 7.7 GeV to 0.15 at 39 GeV. Lower values 
of U2(pt) are observed for anti-protons compared to pro- 
tons at all energies. The difference in the V2 values for 



protons and anti-protons increases with decreasing beam 
energy. The lower panels of Fig. [T] show the pr depen- 
dence of the difference in V2 for protons and anti-protons. 
No significant px dependence is observed, as character- 
ized by the horizontal line fits. The negative values of 
the anti-proton V2 at low pr at ^^snn = 11.5 GeV could 
be due to absorption in the medium . 

The V2{pt) behaviour for A(wds), A(uds) and S~(c?ss), 

S^((iss) is similar to that for protons (uud) and anti- 
protons (uud). In all cases, the baryon anti-particle V2 is 
lower than the corresponding particle V2- The V2{pt) dif- 
ference for A and A is in agreement with previous STAR 
results at ^/snn = 62.4 GeV 0. For the mesons 7r+ (w J) , 
TT~ (ud), and K~^{us), K~(us), the differences are smaller 
than those for the baryons ( the anti-particle convention 
from [2^ is used for mesons). At ^/snn = 7-7 GeV, the 
V2{pt) difference between and is a. factor 5-6 
smaller as compared to the baryons, with having a 
systematically larger V2{pt) than the K~ . On the other 
hand, the V2{pt) of the tt^ is larger than the V2{pt) of 
the TT"*". However, the magnitude of the difference for 
pions as a function of energy is similar to that for the 
kaons. The details of the pr dependence of the difference 
in V2 between particles and corresponding anti-particles 
can be found in Ref. [20j . 

Figure [2] summarizes the variation of the px indepen- 
dent difference in V2 between particles and corresponding 
anti-particles with ^/snn. Here, V2{X) — V2{X) denotes 
the horizontal line fit values of the difference in V2 (pt) be- 
tween particles X (p, A, E!~, tt"*", K'^) and corresponding 

anti-particles X (p. A, S"'^, 7r~, K~). Larger V2 values are 
found for particles than for antiparticles, except for pions 
for which the opposite ordering is observed. A monotonic 
increase of the magnitude of Av2 = V2{X) — V2{X) with 
decreasing beam energy is observed. The data can be 
described by a power-law function. 

While in Au+Au collisions at ^/snn = 200 GeV a sin- 
gle NCQ scaling can be observed for particles and anti- 
particles, the observed difference in V2 at lower beam 
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FIG. 3. (Color online) The upper panels depict the elliptic flow, V2, as a function of reduced transverse mass, {rriT — mo), for 
particles, frames a) and b), and anti-particles, frames c) and d), in 0-80% central Au-|-Au collisions at ^/snn ~ 11-5 and 62.4 
GeV. Simultaneous fits to the mesons except the pions are shown as the dashed lines. The difference of the baryon V2 and the 
meson fits are shown in the lower panels. 



energies demonstrates that this common NCQ scaling 
of particles and anti-particles splits. Such a breaking 
of the NCQ scaling could indicate increased contribu- 
tions from hadronic interactions in the system evolution 
with decreasing beam energy. The energy dependence 
of V2(X) — V2{X) could also be accounted for by con- 
sidering an increase in nuclear stopping power with de- 
creasing if the V2 of transported quarks (quarks 
coming from the incident nucleons) is larger than the 1)2 
of produced quarks [25| . Theoretical calculations [26j 
suggest that the difference between particles and anti- 
particles could be accounted for by mean field potentials 
where the K~ and p feel an attractive force while the K'^ 
and p feel a repulsive force. Most of the current theoreti- 
cal calculations can reproduce the basic pattern, but fail 
to quantitatively reproduce the measured V2 difference 
[2J-[23|- So far, none of the theory calculations describes 
the observed ordering of the particles. Therefore, more 
accurate calculations from theory are needed to distin- 
guish between the different possibilities. Other possible 
interpretations for the observation that the tt^ V2{pt) is 
larger than the 7r+ V2{pt) is the Coulomb repulsion of 7r+ 
by the mid-rapidity net-protons (only at low pr) and the 
chiral niagnetic effect in finite baryon-density matter [28j . 
In Ref. [201 , the study of the centrality dependence of Av2 
for protons and anti-protons is extended to investigate, 
if different production rates for protons and anti-protons 
as a function of centrality could cause the observed dif- 
ferences. It was observed that the differences, Av2, are 
significant at all centralities. 

The V2{mT — mo) and possible NCQ scaling was also 
investigated for particles and anti-particles separately. 
Figure [3] shows V2 as a function of the reduced trans- 
verse mass, (tot — mo), for various particles and anti- 
particles at y/sNN ~ 11-5 and 62.4 GcV. The baryons 
and mesons are clearly separated for ^stvat = 62.4 GeV 



at {rriT ~ rrio) > 1 GeV/c^. While the effect is present 
for particles at = 11.5 GeV, no such separation is 

observed for the anti-particles at this energy in the mea- 
sured (mt — Too) range up to 2 GcV/c^. The lower panels 
of Fig. |3] depict the difference of the baryon V2 relative to 
a fit to the meson data with the pions excluded from 
the fit. The anti-particles at ^/snn = 11.5 GeV show a 
smaller difference compared to the particles. 

In Fig. m the U2(tot — too) values scaled on both axes 
with the number of constituent-quarks are presented for 
^/sNN = 11.5 and 62.4 GcV. A simultaneous fit [2^ to all 
data sets at a given energy is shown as the dashed black 
line. The ratio of the data to these fits is shown in the 
lower panels. Most of the ratio values arc within ±10% of 
unity, showing that the NCQ scaling holds for the group 
of particles while for the group of anti-particles the mea- 
sured kinematic range is not enough to conclude. The 
mesons are an exception. At the highest (tot — TOo)/ng 
values, the meson data point for ^ysNN = 11-5 GeV 
{pt = 1.9 GeV/c) is 2.3a lower than those of the other 
hadrons. This is comparable to the observed deviation 
at y^J^N = 7.7 GeV {pr = 1.7 GeV/c) by I.Sct The 
smaller V2 values of the (/)(ss) meson, which has a smaller 
hadronic interaction cross section [30| , may indicate that 
hadronic interactions become more important than par- 
tonic effects for the systems formed at collision energies 
< 11.5 GeV [311,133. 

In summary, the first observation of a beam-energy de- 
pendent difference in V2{pt) between particles and corre- 
sponding anti-particles for minimum bias ^/s'nn = 7.7- 
62.4 GeV Au+Au collisions at mid-rapidity is reported. 
The difference increases with decreasing beam energy. 
Baryons show a larger difference compared to mesons. 
The relative values of f 2 for charged pions have the oppo- 
site trend to the values of charged kaons. It is concluded 
that, at the lower energies, particles and anti-particles are 
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FIG. 4. (Color online) The number-of-constituent quark scaled elliptic flow {v2 / nq){{mT — mo)/nq) for 0-80% central Au+Au 
collisions at ^/snn = 11.5 and 62.4 GeV for selected particles, frames a) and b), and corresponding anti-particles, frames c) 
and d). The dashed lines are simultaneous fits [2^ to all of the data sets at a given energy. The lower panels depict the ratios 
to the fits, while a ±10% interval is shown as the shaded area to guide the eye. Some data points for (f> and H are out of the 
plot range in the lower panels of frames a) and c). 



no longer consistent with the single NCQ scaling that was 
observed for -y^s^N = 200 GeV. However, for the group 
of particles the NCQ scaling holds within ±10% while for 
the group of anti-particles the difference between baryon 
and meson V2 continues to decrease to lower energies. We 
further observed that the meson V2 at the highest mea- 
sured rriT — rriQ value is low compared to other hadrons 
at y^SNN = 7.7 and 11.5 GeV with l.Sa and 2.3cr, re- 
spectively. 
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